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We study the Kondo-Heisenberg model on the honeycomb lattice at half-filling. Due to the 
vanishing of the density of states at the fermi level, the Kondo insulator disappears at a finite 
Kondo couphng even in the absence of the Heisenberg exchange. We adopt a large-N formulation 
of this model and use the renormalization group machinery to study systematically the second 
order phase transition of the Kondo insulator (KI) to the algebraic spin liquid (ASL). We note that 
neither phase breaks any physical symmetry, so that the transition is not described by the standard 
Ginzburg-Landau- Wilson critical point. We find a stable Lorentz-invariant fixed point that controls 
this second order phase transition. We calculate the exponent u of the diverging length scale near 
the transition. The quasi-particle weight of the conduction electron vanishes at this KI-ASL fixed 
point, indicating non-Fermi liquid behavior. The algebraic decay exponent of the staggered spin 
correlation is calculated at the fixed point and in the ASL phase. We find a jump in this exponent 
at the transition point. 



I. INTRODUCTION 

The interplay between charge and spin degrees of free- 
dom has been a focus of research in complex materials 
such as cuprates and heavy- fermions.^ The clearest ex- 
ample of this interplay is the quantum critical point seen 
in many heavy-fermion materials. At zero temperature, 
the heavy Fermi liquid (HFL) phase disappears exactly 
at a point where the anti-ferromagnetic (AF) magnetic 
ordering grows. Non- Fermi liquid behaviors are seen in 
the quantum critical region, i.e. the V-shaped region 
above this critical point. ^ The theoretical understand- 
ing of why these two seemingly different phases, i.e. AF 
ordered and heavy fermi liquid, should collapse at one 
point and a clear understanding of the non-Fermi liquid 
behaviors in the quantum critical region are poor at the 
moment. 

It is by now understood that the spin density wave 
approach to understand the quantum critical point 
in heavy-fermions, known as the Moriya-Hertz-Millis 
theory^, only accounts for small deviations from Fermi 
liquid theory. New theoretical approaches for under- 
standing the quantum critical heavy-fermions arc needed. 

In an attempt to find an alternative for the Moriya- 
Hertz-Millis theory, Senthil et al. proposed recently that 
the AF-HFL transition might be controlled by an unsta- 
ble spin liquid fixed point.''''* This picture is very similar 
to the Deconfined quantum critical point in the context 
of the second order phase transition between Neel and 
valence bond solid (VBS) ground states.^ In that case 
the transition to VBS happens due to the existence of an 
unstable spin liquid on the magnetically disordered side 
of the critical point. Deconfined quantum critical points 
open up the possibility of second order phase transitions 
between "unrelated" phases like AF and HFL. It also 
has the advantage of giving a plausible scenario for the 
non-Fermi liquid behavior seen in experiments in heavy 
fcrmion materials. This view has had some successes in 
explaining some of the experimental observations.^''' In 
this paper, we further explore this point of view. 



In search for a microscopic model that provides this 
type of quantum critical point, we study the Kondo- 
Heisenberg model on the honeycomb lattice at half-filling. 
The heavy fcrmion quantum critical point in this model is 
simplified. It corresponds to a point, where both charge 
gap and spin gap vanishes. This model is given by the 
following Hamiltonian {Jk > 0, Jh > 0): 



H 



-t 



E 



H.A +Jk J2 S^■S^ + JH ^^'^J 



(1) 

where i and j live on the honeycomb lattice sites and 
a is the spin index: {t,i}- The S; and Si denote the 
conduction electron spin and the localized spin at the 
site i respectively. What makes the honeycomb lattice 
very interesting to study is the fact that the Kondo gap 
vanishes at a finite coupling constant even as Jh — *■ 0. 
This is due to the fact that, in contrast to the square 
lattice, the density of states of the conduction electrons 
vanishes at the Fermi level.® 

The Kondo-Heisenberg model on the honeycomb lat- 
tice at half-filling is also interesting, because it can be 
studied by quantum Monte Carlo without the fcrmion 
sign problem. ^'"'^'^ The present paper can be considered 
a prelude to such a study. In this connection wc men- 
tion the quantum Monte Carlo study of the Kondo lat- 
tice model {Jh = 0) on the square lattice by Assaad.*° 
He found a quantum phase transition of the AF ordered 
ground state to a disordered one caused by the Kondo 
exchange. The quasi-particle gap however did not van- 
ish at the transition point, i.e. the magnetic transition 
happened inside the KI phase. This made the transi- 
tion a magnetic transition belonging to the 0(3) nonlin- 
ear sigma model universality class. Note that the nested 
Fermi line of the tight-binding band provides the insta- 
bility towards AF at (tt, it) even in the absence of the 
Heisenberg exchange Jh ■ In contrast the honeycomb lat- 
tice tight binding band has Fermi points in the corner of 
the Brillouin zone called the Dirac points. This causes 
the nesting to be extremely weak as compared to the per- 
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FIG. 1: A scenario for the zero temperature phase diagram 
of the Kondo-Heisenberg Hamiltonian on the honeycomb lat- 
tice at half-fiUing [given by Hamiltonian Eq. (1)] as a function 
of Jk (t = 1) for a sufficiently small Jh- In contrast to the 
same model on a square lattice, Kondo gap denoted schemat- 
ically by the dotted-dashed line vanishes at a finite coupling 
constant J^k- The weak nesting of the Dirac nodes of the 
honeycomb lattice makes the KI transition to a spin liquid 
ground state plausible. 



fectly nested lines of the Fermi line in the square lattice. 

Next we consider different scenarios for the ground 
state phase diagram of our model, where we consider 
varying Jk {t = 1), keeping Jh to be small. In the 
limit Jk <JC Jh, the Heiscnberg exchange dominates and 
will lead to Neel ordering with opposite spins on the AB 
sublattices of the honeycomb lattice. 

In the first scenario, sketched schematically in Fig. 1, 
the finite coupling Kondo transition is to a magnetically 
disordered ground state (spin liquid). In this scenario 
the AF ordered ground state docs not survive up to the 
KI transition point for sufficiently small Jh ■ The Kondo 
physics has stabilized a spin liquid ground state in the 
region between the KI and AF ground states. The charge 
gap closes at Ji^, as the KI gives way to a semi- metal. 
The spin gap may or may not be finite, depending on 
the nature of the spin liquid. As Jk is further reduced 
a transition to the Neel phase occurs. This scenario will 
of course be very exciting, if a spin liquid can indeed be 
stabilized as the ground state over some parameter range. 

The other exciting possibility is that a continuous tran- 
sition between two distinct orders (AF and KI) exists 
over a finite parameter range. This scenario is sketched 
schematically in Fig. 2. In this scenario SL is unstable 
to Neel ordering. This will be an example of the decon- 
fined quantum critical point mentioned earlier. Finally, 
scenarios in which an overlapping region, where a Neel 
state coexists with the Kondo insulator (similar to the 
square lattice result of Assaad^°); or a first order tran- 
sition can not be ruled out. With these motivations we 
now sketch an outline of this paper. 

In this paper we begin by a mean-field study of the 
Kondo-Heisenberg model Eq. (1) on the honeycomb lat- 
tice. We adopt a fermionic representation for localized 
spins. We make a mean-field decoupling of Si ■ Sj in 
the fermionic hopping channel. This is sometimes re- 
ferred to as the resonating valence bond (RVB) decou- 
pling. The Kondo interaction is decoupled in the usual 
mean-field way with the hybridization "order parameter" 
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FIG. 2: Another scenario for the zero-tempreature phase dia- 
gram for the Kondo-Heisenberg Hamiltonian at half-filling on 
the honeycomb lattice for some (maybe "fine tuned" ) Heisen- 
berg exchange Jh as a function of Jk- In contrast to the 
same model on a square lattice, both the quasi-particle gap 
denoted schematically by the dotted-dashed line and Neel or- 
der paramater denoted by the solid line vanishes at Jk = Jk ■ 



b = -fi^ '^ij ■ "^^^ nonzero b will correspond to the 
KI phase, where the gap in the dispersion is proportional 
to b. For any Jh > we find a continuous transition be- 
tween the KI (6 ^ 0) and a spin liquid state (see Fig. 3). 
This spin liquid is characterized by Dirac fermions cou- 
pled to U(l) gauge field, and has been called the algebraic 
spin liquid (ASL), because the spin correlation decays as 
a power law and spin excitations arc gapless.^^ 

As mentioned earlier, in the Jk ^ Jh hmit, localized 
spins will be Neel ordered. Our mean field decoupling can 
not capture this.^^ Here wc proceed with the assumption 
that the phase diagram is given by Fig. 1 and our main 
interest is to study the KI to SL transition. For this pur- 
pose the RVB decoupling is a reasonable starting point. 




FIG. 3: The mean-field phase diagram of the Hamiltonian 
given by Eq. (1) in the Jk — Jh space {t = 1). For any Jh > 0, 
we find a transition from the KI to the algebraic spin liquid. 

The rest of the paper is devoted to studying the critical 
properties of the fixed point that controls the transition 
from the KI to ASL. Wc note that neither phase breaks 
any physical symmetry, so that the transition is not de- 
scribed by the standard Ginzburg-Landau- Wilson critical 
point. What changes at the transition is the dynamics 
of an emergent gauge field. The emergent gauge field is 
confined in the KI phase and it is deconfined in the ASL 
phase. 

In Sec. II B wc generalize our model by letting the spin 
indices to run from 1 to N rather than just "up" and 
"down" . The saddle point approximation (i.e. the mean- 
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field) becomes exact as iV ^ cxj. In Sec. Ill we derive 
tlie low-energy Lagrangian density which sets the stage 
for a systematic expansion. 

Section IV develops the propagator for the Kondo field 
b, and the gauge field in the leading order. The Kondo 
field propagator can be tuned to a masslcss point at Jk = 
J'^. In the J > regime the b field condenses to \b\ 
and we will be in the Higgs phase where the gauge field 
is massive. This is nothing but the Kondo physics, where 
the excitation gap is proportional to In the J < 
regime the Kondo field is not condensed, and is massive. 
The physics is described by Dirac fermions coupled to 
a U(l) gauge field. This field theory, known as QED3, 
has been studied extensively and is understood to 
be an algebraic spin liquid due to algebraic correlation of 
spins. This phase is believed to be deconfined.^'^ 

Section V A contains the main technical calculation of 
this paper. Wc calculate the exponent which describes 
the diverging length scale of the transition. Wc consider a 
relativistic fixed point, which allows us to adopt standard 
field theory methods. In Sec. VB we calculate the de- 
cay exponent of the staggered localized spin. In Sec. VI 
we show that the relativistic fixed point is stable and 
is therefore the appropriate fixed point to study for the 
transition. 



II. FINITE COUPLING KONDO TRANSITION 



Mean-Field Transition 



together with 



(6) 



will result in the SU(2) commutation relations for spin 
S. However, the Hilbcrt space of the localized spin Si is 
restricted to only two elements: {\iti) , l-IJ-i)}- Therefore 
identifying localized spin by fermionic operators given by 
Eq. (6) has to be accompanied by the constraint 



1. 



(7) 



Combining Eqs. (3) — (7); together with the following 
identity for the Pauli matrices 



fab • fed — 25adSbc — Sabred 



(8) 



one finds a fermionic representation of the spin-spin in- 
teractions in the Hamiltonian of Eq. (2): 



Si ■ Sq 



H.c. 
H.c 



1 

4' 
1 

^4' 



(9) 
(10) 



where the sum over the spin indices is understood. 

The four-fermion interaction terms of Eqs. 9 and 10 
will be replaced in the mean-field (MF) simplification by 
a quadratic interaction: 

A = A,A2 ^ (Ai) A2 + ii (A2) - (A,) {A2) , (11) 



Our starting point is the Hamiltonian 

= ^ [cfc'^ + H.c}j +Jk s^■S,+JH S,-S,, 

{ij),a i (ij) 

(2) 

where i and j live on the honeycomb lattice sites and a 
is the spin index: {t,i}- The Si denotes the conduction 
electron spin at the site i and the capital Si denotes the 
localized spin at the site i and they both are SU(2) spins 
with an SU(2) spin algebra. The conduction electron spin 
is given by 



1 



at b 



(3) 



where A denotes a generic four-fermion interaction term. 
To make the interaction quadratic one has to ignore (^i — 
{Ai)){A2 — (^4.2)). So the mean-field parameters have to 
satisfy the self-consistency condition: 



dF, 



MF 



d{A,) 



0. 



(12) 



We make the interaction terms given by Eqs. (9) 
and (10) quadratic by choosing the following mean-field 
parameters: 



(13) 



where cr = (cTj;, Uj,, ct^) arc the Pauli matrices. Wc adopt 
a fermionic representation for localized spins, where their 
Hilbert space {|f|"i) , l^i)}, is constructed using Fermionic 
operators: 



l^.)=//^|0), 
l^.)=.//^|0). 

The anti-commutation relation: 



(4) 



(5) 



\Xij\ 



(14) 



The first parameter |6i|e*^' lives on sites. The second pa- 
rameter |xij|e*°'^ lives on links. Wc assume that these 
parameters have a constant magnitude throughout the 
lattice. When they take a non-zero value, small devia- 
tions from their "constant" magnitudes cost energy: they 
are locally stable. \bi \ = b ^ Q corresponds to the Kondo 
insulator phase, where the gap in the excitations is pro- 
portional to b. Wc also note that the conventional de- 
coupling of Si ■ Sj leads to the anti-ferromagnetic order 



parameter (Siz). In this paper we adopt the alternative 
decoupling of Eq. (14), which is often called the RVB 
(Resonating Valence Bond) decoupling. RVB decoupling 
leads to a spin liquid state. 

The U(l) gauge freedom in defining / gives the freedom 
to eliminate 6i by the following gauge transformation: 



/r 

a, 



(15) 



This transformation leaves the total flux modulo 2tt 



aij (mod 27r) 

O 



(16) 



through any closed loop invariant. Next we make a choice 
Xij = X where x is real. This choice corresponds to the 
case where the total flux through each hexagon modulo 
27r is zero. 

Our mean-field choices result in the following mean- 
field Hamiltonian: 



JhX 



H.c. 



(17) 



^f[JKh'+'^^x' 



where M is the number of unit cells. The different co- 
efficients of the constant terms are due to the fact that 
there are 2 sites as opposed to 3 links per unit cell. We 
also mention that, at the mean- field level, the constraint 
given by Eq. (7) is enforced on average by having a chem- 
ical potential for / fermions. This chemical potential is 
zero due to f-fermions particle-hole symmetry. 

From Hm f , the MF free energy can be calculated nu- 
merically. MF parameters satisfying the self consistency 
condition are then obtained. The numerics confirm the 
finite coupling Kondo transition in a background Xi for 
any Jh > 0. Next we show this finite coupling Kondo 
transition analytically and find the exponent (3 



\b\ cx (J 



K 



jcr\P 



(18) 



by focusing on the low-energy theory near the Dirac 
nodes. 

We focus first on the c electrons kinetic term 
—tY^ (c^^Cj + H.c^ . We do the Fourier transformation 



CA{k) = j^Y.^""''^^(^) 



(19) 



R 



where R are the Bravais lattice vectors and ca{R) is 
the conduction electron annihilation operator at the A 
sublattice site of the unit cell positioned at R. The sim- 
ilar transformation is done to get CB{k). Conduction 





FIG. 4: The Brillouin zone of the honeycomb lattice. Corners 
of the Brillouin zone are where the tight-binding gap vanishes. 
The independent low-energy modes are denoted by the fiUed- 
circles around the two Dirac nodes ztfeo- 



electrons kinetic term in momentum space will take the 
form: 

-^telcj+H.c. = -tJ2^ik)c\{k)cB{k) + H.c., (20) 



where 



r(fc) = 1 + e''"''-'^^"^^'' + (2»7i+»7: 



'2) 



(21) 



and rji ~ ^(1,0) and r]2 = /(— 1/2, ■\/3/2) are two of the 
three vectors that connect A sublattice sites to nearest 
neighbor B sublattice sites. / is the length between near- 
est neighbor sites which will be set to the unit length: 
I = 1. 

T{k) vanishes at the Dirac nodes ifcu, which is given 
by 



k, 



(0, 



An X 



(22) 



Near the nodes as schematically represented by the 
filled circles in the Fig. 4, r(±A;£) + q) is given by 



(23) 



To write the conduction electron kinetic term near the 
Dirac nodes, we use the following notation, which will 
also be helpful in Sec. Ill: 



V'±(g) 



CA{±kD + q) 
CB{±kD + q) 



(24) 



After collecting Eqs. (20) — (24), we get 



^tcjcj + H.c. - — ^7/'+(9)t(qicr2 + 92cri)V'+(g) 
(ij) 1 

M 



X^V-'- (9)^(91^^2 - q2<Ji)tp-iq): 

(25) 



where "r^" is used, since in contrast to Eq. (20) the mo- 
mentum sum is restricted to be near the Dirac nodes. 
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Next we obtain the energy dispersion near the node 
The MF Hamiltonian in momentum space and 
near the node denoted by H+{q), is given by 



Mq) 



(26) 



where we have dropped the constant term of Eq. (17) and 
we have used a similar notation as in Eq (24) for the / 
fcrmions: 



-Xq) 



fA{±kD+q) 
fB{±kD + q) 



The 4x4 matrix 7i+ is given in block form by 

n+ = 



2vc{qia2 + q2Cfi) 
JKb 



2vf{qia2 +920-1) 



where Vc and Vf are given by: 



3 

-t, 



(27) 

(28) 

(29) 
(30) 



It will be clear why notations Vc and v/ is used shortly. It 
is straightforward to diagonalize the above matrix. The 
four eigenvalues are given by: 



e{q) = ± ^{v, - Vf)\q\ ±^{vf+ v,)W + Jlb^ j /2. 

(31) 

We have dropped the subscript + from e{q), since the 
dispersion near the other node is the same. 

To understand the above dispersion better, we mention 
its behavior in different regimes: 

• 6 = 0. The MF Hamiltonian has two decoupled 
parts: kinetic terms for the conduction electrons 
and for f fermions. In this case the 4 eigenvalues 
near the nodes become ±Wc|<z| and ±u^|g|. They 
are the massless Dirac dispersions for conduction 
electrons c and fermions / with velocities Vc and Vf 
respectively. The phase with 6 = is an algebraic 
spin liquid phase as will be clear in Sec. VB. This 
is closely related to the work of Rantner and Wen, 
where Dirac fermions coupled to a gauge field was 
found to be an algebraic spin liquid. 

• Vf ~ Vc,b ^ 0. The dispersion is the massive rela- 
tivistic Dirac dispersion 



eiq) 



±\/v1q^ 



{JKh/2f 



(32) 



where each eigenvalue is doubly degenerate. JKb/2 
is the mass of our relativistic quasi-particles. In 
Sec. Ill, we will see that the low-energy theory, in 
the limit Vf = Vc, is Lorentz invaraint. The above 
relativistic dispersion is a sign of this Lorentz in- 
variancc. 



sign of Vf . The sign of x (f^f) has dramatic con- 
sequences on the shape of the dispersion and also 
on the continuum formulation of the model near 
the phase transition. Due to the coupling between 
f-fermions and the conduction electrons the free en- 
ergy is not invariant under the mapping x i— > — X- 
In the Kondo insulating phase (6 7^ 0); the disper- 
sion for Vf > is gapped. However the dispersion 
for Vf < has a gaplcss ring around each Dirac 
node with the radius 



JKb 



2y/-VcVf 



(33) 



It turns out that the positive (self consistent) x 
{vf > 0) is more stable. This is not surprising from 
the above fact that the dispersion for Vf > is 
gapped as opposed to Vf < which is gapless. We 
arc only concerned with the stable solution x > 
in this paper. 

The MF transition at zero temperature can be ana- 
lyzed by minimizing the total energy 



E^4J2 <q)+^(JKb' + ^fx' 



e(q)<0 



dE 
'db 



0, and 



9E 
dx 



(34) 



(35) 



The coefRecient 4 = 2 x 2 is for counting spins and the 
two nodes. The condition ^ = at 6 = results in the 
self consistent Xi which varies slowly across the transi- 
tion. For Vf > 0, the negative energy bands arc given by 

(±{vc - Vf)\q\ - v/(w/+Wc)2g2 + Ji6^) /2. Therefore 



s(q)<0 



(36) 



Combine Eqs. 34 — 36, to arrive at the self-consistency 
equation for b: 



Jk 



^^vw+vcrq' 



Jib- 



1. 



(37) 



It is clear from the above equation that the Kondo phase 
does not survive down to Jk = 0, in contrast to the case 
when the conduction band possess a Fermi surface. The 
Kondo insulator phase disappears at a finite coupling 
constant JJ^. Taking the limit A/" 00, we find the 
following dependence of the Kondo parameter \b\ as a 
function of the distance from the transition Jk — Jk'- 



j: 



\b\ (X (Jk - J'kI, 

Vc + Vf 



(38) 



K 



A 



where A ^ 1 is a large-momentum cutoff around the 
Dirac nodes ±fc^. 

To summarize, the MF solution has two interesting 
features: 
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• The Kondo insulator phase disappears at a finite 
couphng constant. 

• The Kondo condensate \b\ vanishes hnearly as the 
distance from the critical point. We emphasize that 
a Landau-type expansion for the b mean-field pa- 
rameter would yield the result |6| oc {Jk — J'j^Y^^ 
near the transition point. The appearance of non- 
trivial critical exponent at the mean-field level 
is characteristic of the coupling of b to gapless 
fermions which is addressed systematically in the 

expansion. 



B. Large N Formulation : Mean-Field Justified 

In this section we extend the physical model with two 
spin-flavors to a generalized model with N spin-flavors.^^ 
We show that the mean-field solution is the stable solu- 
tion in the infinite N limit. 

The conduction electrons kinetic term 
^ ^c°^c° + iJ.c.^ is already in the generalized form. 

We just let the spin index a to run from 1 to N. The 
spin-spin interactions as written in Eqs. 9 and 10 is 
also easily generalized by letting the spin indices to run 
from 1 to N. In addition we control these interactions 
by replacing i with . So the generalized version of the 
Hamiltonian of Eq. (2) is given by 



H = -t 

_ Jk_ 
2N 



H.c. 



2N 



Y.[icffnif-'ct) + H.c 

i 



(39) 



where the sum over spin indices is understood. The above 
Hamiltonian has to accompanied by the local constraint: 



/ J J i J i try ' 



(40) 



The above generalized model will coincide with the 
Hamiltonian Eq. (2) for N=2. In fermion coherent 
state representation, using Grassmann variables, the Eu- 
clidean Lagrangian and the partition function take the 
form 

Li,([c,/,ao]) = ^^c^(xo)aoc^(xo) + i?(a:o)9o/r(a:o) 
-ia.Q /r(xo)/f(xo) - - 




Z = Vlccffao] e' 



' f dxoLEixo) 



where Xq denotes the imaginary time. The field Oio, when 
integrated over, will produce delta functions at each site 
enforcing the local constraint given by Eq. 40. The no- 
tation ttio is used because it will serve as the time com- 
ponent of the emergent gauge field. 

Next we do the standard Hubbard-Stratonovich trans- 
formation to decouple the four-fermion terms in the ac- 
tion at a cost of introducing the complex fields bi, which 
live on sites and Xiji which live on links: 

Le{[c, /,6,x,ao]) 
= Y.c'^doc'^-tY,{c-c^ + c.c-) 

* (ij) 



J2{fndo-ia,o)ft+iNa,o/2 



(42) 



Z = I 'D[ccffb*bx*xao]e' 



The xo dependence of Grassmann and complex fields in 
the Lagrangian is understood as well as the dropped in- 
dices in the measure of the path integral. 

The action is quadratic in the Fermionic fields and it 
can be integrated out at the cost of obtaining an effective 
action for bi and Xij fields, which is highly nonlocal. Since 
there is no mixing of different spin-flavors, the resulting 
effective action will be proportional to N: 



LE{[b,x]) 
N 



= logdctX + i-^ |6,|2 + I^^^l' 



(41) 



(43) 

is a 4// X 4J\f Hermitian matrix, which encodes the 
quadratic interaction of the fermionic fields as given by 
Eq. (42), e.g. Mc,f, = b,.^'^ The "detAl" makes the 
effective Lagrangian in terms of bi and Xij fields nonlocal. 

The saddle point approximation for the above La- 
grangian LE{[b,x]) is exact as iV oo. The saddle 
point equations will be the mean-field equations given by 
Eq. (12), and this is how the mean-field developed in the 
previous section is justified. Furthermore the Lagrangian 
given in Eq. (42) will set the stage for a systematic 1/A^ 
expansion around the mean-field solution. 



III. BEYOND MEAN-FIELD: LOW-ENERGY 
EFFECTIVE FIELD THEORY 

In this section wc obtain the Lagrangian density near 
the critical point starting from the microscopic La- 
grangian. There are two independent nodes ifcjj and 
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near each node there are fermions from A and B sub- 
lattices. We package these c and / fermions into four 
component spinors ^' and <&. This is done in a way to 
make the kinetic terms for both fields look the same. 

Since there is no mixing of different spin-flavors in the 
microscopic Lagrangian Eq. (42), we will drop this index 
throughout the derivation of the Lagrangian density. The 
spin-flavor index will be resurrected at the end. 

The kinetic term for the conduction electrons was dis- 
cussed in Sec. II A: 



c\cj -f H.c. - Vc^'^:+{q)\qi(J2 + q2(Ti)i'+{q) 



q 

E 
1 



(44) 



Note that under the transformation 'ip-{q) (T2ip-{q) 
the second term has the same form as the first one. We 
combine the two 2-component spinors, associated with 
the two Dirac nodes, to form a 4-componcnt spinor ^'(t/): 



*(g) 



■ i'+iq) ' 

cr2'ip-{q) 



(45) 



To write the kinetic term in terms of "^{q), we define the 
4x4 block diagonal matrices r^: 



cr^ 
a. 



I-2X2- 



(46) 



The kinetic term then takes the form 



clci + H.c. 



. J2'^il)Hqn-2 + q2Ti)^{q)- (47) 



After taking the continuum limit, the Lagrangian den- 
sity corresponding to this Hamiltonian is: 



(48) 

To use the field theory tcchniquGS one would write 
in a different form 



where is defined to be 



T3- 



T2d2)^, (49) 



(50) 



Next we focus on {x^oI^ f! + Xljlj ff) term in 
Eq. (42). Consider the mean-field case Xij = X > 0- The 
phase fluctuations will be included as gauge field fluc- 
tuations later on. Ignoring the phase fluctuations, this 
term looks exactly like the kinetic term of the conduction 
electrons given by Eq. (44) , but with an opposite sign: 



H.c. ^ Vf^(j>+{q)^{-qi<T2 - q2cri)(l>+{q) 



^ 4>-{q)\~qia2 + q2<Ji)(f>^{q). 

(51) 



1 



We make this look like Eq. (47) by a further transfor- 
mation to absorb the minus sign and we again combine 
the 2-component spinors near the two nodes to form a 
4-component spinor $(</): 



o-3</>-H(q) 

-icri</)_(q) 



(52) 



On the grounds of gauge invariancc, we write the La- 
grangian density corresponding to the "kinetic" term of 
the / fermions 

Z:^' = l>r3(ao-mo)<I>+w/$(-ri(ai-mi)+T2(52-m2))$, 

(53) 

where the non-compact gauge field a(x) appears as a 
connection to produce covaraint derivative. It is related 
to the the phase aij through the relation a{x) = aijfiij, 
where fjij is the unit vector connecting the site i to the 
site j. We have used the same definition as in Eq. (50): 



(54) 



En route to the continuum limit, we have implicitly 
dropped the compact character of the gauge field. This 
is only valid if the gauge fluctuations are not very strong. 
In the imaginary time evolution, there will be events on 
the lattice scale that change the flux of the gauge field by 
2tt. If these events do not proliferate, we will be in small 
gauge coupling constant regime and taking the contin- 
uum limit is justified. We return to this point at the 
end of this section. 

Now we focus on the interaction term J^i^i'^lfi- We 
write this in the momentum space and restrict c and / 
fermions momentum to be near ikjj, and b momentum 
to be near fc = 0. Simple manipulations 

^liq)MQ') = i^l{q)a3 iasMl')) . (55) 

^Uq)c^-{q') = (a2V-(g))ta3(-*ai0_(q')), (56) 

combined with the definitions for ^' and <f> fields, result 
in 



V'l (q)0+ (g') + (?)(/)_ (?') = vl/t (q)r3<i>(q') 



(57) 



= 'fiqMq')- 

Therefore in the continuum limit the Lagrangian density 
corresponding to J2i (&j4/» + Kfl^ij wiU be 



(58) 



Now collect jC^"*, jC^"* and to write the La- 

grangian density 



5^'''$'^ + + — l&p + • • • , 

Jk 



(59) 



where the spin indices are restored. 7^ is defined to be 

1^= {T3,-gTi,gT2) , (60) 
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and 7p is a shorthand for 7^ with g = 1. The velocity Vf 
in Eq. (53) is absorbed by scaUng x and a. As a result, 
the velocities enter only as the dimensionless ratio Vc/vf 

Vc/Vf 

Via 7p . 

Taking the continuum limit, which is essentially 
a coarse graining process, generates new interactions 
among coarse grained fields.^" The ellipsis in Eq. (59) de- 
notes such new interactions. Not all possible terms are al- 
lowed though. Since the coarse graining process does not 
break symmetries of the microscopic Lagrangian, the gen- 
erated interactions should respect these symmetries. 

Next we set Vf = Vc to make the low-energy theory 
Lorcntz symmetric. The action in the presence of two 
different velocities will not respect the Lorcntz symmetry. 
Under RG transformation — in the absence of Lorentz 
symmetry — the velocities will change, since there is no 
symmetry to protect them from changing. We shall show 
later that the Lorentz invaraince is protected in the RG 
sense, in that small deviation from Vf = Vc will scale 
to zero under RG transformation. Thus we focus our 
discussion on this Lorentz invariant fixed point. 

The Lorentz-symmctric Lagrangian density is given by: 



N .-, 



Jk' 



(61) 



N 



+ ^ (e^^Aa^flA) + Ng\ (a^ + m^) b\^ + 



where is a shorthand for with g = I. The last 
two terms are two examples of the terms that will be gen- 
erated in the coarse graining process. They are written 
on the grounds of gauge invariance. In Sec. VI we show 
that these new terms and other terms grouped in the 
ellipsis are irrelevant at our finite coupling Kondo tran- 
sition fixed point. Embedded in Eq. (61), is the QED3 
Lagrangian desnity 



N 2 



(62) 



QED3 has been studied before as the low-energy theory 
of the staggered flux phase^^"^''. We have used the recent 
rcsult^^ about the irrelavance of instanton events near the 
QED3 fixed point and took the gauge field to be non- 
compact. In our formulation the physical SU(2) theory 
corresponds to N=2 /owr-component spinors. It could 
have been formulated as Nf = 4 iwo-component spinors 
near each node in the Brillouin zone. To go back and 
forth between these formulations one just needs to use 
Nf = 27V. 



IV. 1/N SYSTEMATIC EXPANSION: FIELD 
PROPAGATORS 



The Lagrangian density Eq. (61), derived in the previ- 
ous section, allows us to develop the machinery of 1/iV 
expansion. The propagators for the fields b and are of 



order 1/N in the large N limit. This will alow us to cal- 
culate different propagators and scaling dimensions order 
by order in 1/iV. 

We start with our notations for propagators of the 
fermion fields, ^' and $. We use a single line notation for 
the $ field propagator and a double-line for 'J propaga- 
tor: 



{^^^)iq) = i- 
where Feynman's slash notation 

is used. 



S^'ii^);' (64) 



(65) 



The Gauge Field a^. Our Lagrangian density Ce has 
a QED3 piece 



(66) 



which has been studied before. ^^"^^ We refer the reader 
to the Appendix B of Ref. 14 for the calculation of the 
QED3 gauge propagator. The QED3 gauge propagator, 
in the leading order, is given by: 



{af,a^) (q) = ^i- 



N\q\ 



(e-i) 



(67) 

where ^ is the Fadeev-Popov gauge fixing parameter. 
^ = 1 corresponds to the Feynman gauge and ^ = 
to the Landau gauge. The \q\ dependence comes from 
calculating the self energy 



o 



(68) 



The charge e has disappeared form the final result. 
This charge appears if one includes the q^ dependence 
in the denominator of the gauge propagator. But in 
the small q (long distance) limit, the \q\ dependence 
dominates over q^ term, no matter how large e is. 
Another way to say this is that the QED3 is self-critical. 
It will flow to a stable fixed point — even if we start with 
a large gauge charge e — at least for some large TV. 

The gauge propagator of our theory, in the leading 
order, will come solely from the f-fermion bubble, given 
diagrammatically in Eq. (68). The b propagation does 
not contribute to the gauge field's self energy in the 
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leading order. Therefore, in the leading order, {af^ai,) (q) 
of our Lagrangian density Eq. (61) is the same as QED3. 

The Kondo Field b. The b field propagator, in the 
leading order, is given by: 



k + q 



{bb*)-' {q) ^ G-\q) = G,, 



(69) 



where G^^iq) = N/Jk + 0(<?^). The self energy 
is easily obtained to be: 
k + q 

, r d^k 

k 



N. , 



2NA 



Tr 



This will result in the following b propagator: 

g 4 



(bb*) {q) 



N(\q\+mb) ' 



(70) 



(71) 



where mt is the effective mass of the 6 propagator, and 
is given by: 



■nib 



4 



'4 « (Jk 



Jk), 



(72) 



is of course 



where J|r is given by = 7r^/(2A). 
nonuniversal, i.e. it depends on how we regularize the 
integral in Eq. (70). However the final result that 



nib oc (Jj^ - Jk) 



(73) 



is scheme independent. 

The above result, that the mass of the b field can be 
tuned to zero, serves as a sanity check for the mean-field 
result of the finite coupling Kondo transition obtained in 
Sec. II A. 

Since the self energy has \q\ dependence, in the long 
distance (small \q\) limit, it dominates over the q^ de- 
pendence coming from The q^ dependence is ig- 
nored in the final expression given in Eq. (71). In the RG 
language, the term is irrelevant at the fixed point. 

We also mention a further notation for the massless 
b propagator. We have used the "double-dashed" hue 
notation for b propagator away from the critical point. 
We use a "simple-dashed" line for the b propagator at 
the critical point: 



(74) 



N\q\ 



Except for the calculation of 7;,, done in Sec. V A 1, mass- 
less b propagator of Eq. (74) is used throughout the pa- 
per. 



The b propagator given by Eqs. (71) and (74) is gauge 
invariant. This will not be the case to all orders in pertur- 
bation theory. The reason is as follows. There is a gauge 
freedom in defining Given that the action is gauge 
invariant, the interaction 6^'^$'^ -|- 6*|>a\|/'» dictates that 
gauge transformation 



$(a;) -> e'"("^)$(a;) 



has to be accompanied by 

b{x) e 



-ia{x) 



b{x), 



(75) 



(76) 



since ^ can not tolerate any gauge freedom. So the b 
propagator in principle should have a gauge dependence. 
However to leading order, the gauge field propagator does 
not enter in &'s self energy and that is why the leading 
order b propagator given by Eq. (71) is gauge invariant. 



CRITICAL PROPERTIES OF THE FIXED 
POINT 



Divergence of The Correlation Length: 
Exponent i' 



the 



In this section we go through the calculation of the 
divergence of the correlation length ^ near the fixed point 
by calculating the exponent v. 



1 



Jk - J] 



(77) 



if I 



We mention first what the correlation length means by 
looking at the Kondo-Heisenberg Hamiltonian in differ- 
ent regimes. At Jk = the Kondo-Heisenberg Hamil- 
tonian has two decoupled parts: the kinetic term for the 
conduction electrons and the anti-ferromagnetic Heisen- 
berg exchange of the localized spins. There is no en- 
tanglement between conduction electrons and the local- 
ized spins in the ground state. In the large Jk limit 
the Kondo term dominates over the Heisenberg exchange. 
The ground state is the product of conduction electrons 
and the localized spin singlets at each site. Since we are 
at half filling every conduction electron is taken away by 
a localized spin. There is entanglement in the ground 
state. But the entanglement is restricted to one site. As 
we approach the critical point, the "size" of these Kondo 
singlets grows and it eventually diverges at the critical 
point. The scale associated with these Kondo singlets is 
what we mean by the correlation length ^. 

The relevant flow of the mass term |6p will introduce a 
diverging length scale ^b near the transition point. Since 
the Kondo phase is charactarized by the condensation of 
the b field, the correlation length ^ is proportional to ^6. 
This results in 



1/ = Vb. 



(78) 
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We obtain vi, by using the scaling relation 



20 



where ^^^(O, mti) is given diagramatically by 



lb 



2-Vb 



(79) 



where 7t, and 2 — rjt characterize the power-law behavior 
of two-point function Gb{k,mh) 



{h{x)b*{x')) = 



(Pk 
(2^ 



^zfc.(.-.')G,(fc,mfc) 



(80) 



S,i/^(0,m6) = 




(84) 



in two different limits: 



06(0,7716 0) cx m^^^\ 
G6(/s^0,0)cx |fc|-2+'7i.. 



(81) 
(82) 




^— 15^7776 log(777h//i). 



It should be noted that since b is not gauge invariant, 7^ 
and 77ft are not gauge invariant. However v by definition, 
given by Eq. (77), is a gauge invariant quantity. Since it 
is not a priori obvious that the scaling relation given by 
Eq. (79) is gauge invariant, we calculate 7^ and rjb in a 
general gauge, and show that v derived from Eq. (79) is 
indeed gauge invar aint. 

The methods we use to calculate 7^ and rjb are very 
different. The calculation of 7^ is done by directly calcu- 
lating the Green's function 06(0,7776). Calculating 77b, in 
this manner, is much more involved. Instead, we calcu- 
late rjb by relating rjb to the scaling dimension of ^"(J)". 
The scaling dimension of is then obtained by per- 

turbing the Lagrangian density Ce Ce + (5^°$° and 
applying the Callan-Symanzik equation for the propaga- 
tor {^^f^^"), in the presence of the perturbation, to ob- 
tain the flow of the composite operator ^°<I>°, and thus 
extracting its scaling dimension. 



This results in 



1- 76 



G^\Q, mb) = (1 - log(7^6//i) 



OC 777. 



.114-1) 



(85) 



from which we obtain 



(86) 



The integrals are evaluated in the dimensional regu- 
larization scheme. ^'^■^^ The mass scale jj, in this scheme 
is introduced to keep track of dimensions in d dimen- 
sional integration and plays the role of a UV cutoff. A 
simple pole in the dimension d reflects the logarithmic 
divergence of the original 3 dimensional integration. We 
only keep track of the IR divergent parts. As far as the 
universal properties of the fixed point are concerned this 
regularization would give the same answer as a cutoff 
regularization. 

We briefly discuss the diagrams involved in 
^^''^(0 , 7776). The first diagram of Eq. (84) has no 
rub dependence. The second and third diagrams are 
equal and are given by: 



To calculate 7f, we calculate Gb(0, 777^), i.e. the uniform 
part of the b propagator slightly away from the critical 
point. From the calculations in the previous section lead- 
ing to Eq. (71) we have 




p + q 



-^rnb\og{mb/ji). (87) 



G^\0,mb)^N 



nib 



1:^^(0,7776), 



(83) 

The last diagram of Eq. (84) contains the vertex that 
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couples b bosons with the gauge field: 




invariance. For the propagator under consideration it 
will read 



(*'^(e^a;)$'^(e^a;')) - 0, 



d . dCra d 

(94) 

where A$ and Aij, are scaling dimensions of the fields <& 
and ^. The set C includes the coupling constants that 
flow. For the perturbation given by Eq. (93) it has only 
1 element:^® 



In the calculation of £^^(0, m) the momentum entering 
this vertex is zero. The result of this integral for fc = is 



Tr 



-9m 



Having the result of the above integral, we get 

q 



(89) 



pi + q 



P2 + q 



E 



pi + q 




= 4jTOfclog(mb//i). 



P2 + q 



(90) 



2. r,b 



C^{5}. 



(95) 



It is more convenient to apply the Callan-Symanzik 
equation to the momentum-space Green's function 
(4'''<I>'') (fc) 



(Pk 
(2^ 



e'''^''-'''^ {■^''{x)^''{x')) . (96) 



Eq. (94) results in 

•I + (3 - A»- A,, + 

(97) 

where "3" is coming from the scaling of the measure cPk 
in the integrand of Eq. (96). 

To apply the Callan-Symanzik equation we first need 
to obtain A$ and A,p. The gauge dependant A$ is ob- 
tained by calculating the propagator 



?76 is easily related to the scaling dimension of the b 
field at the critical point [See Eq. (82)] : 

7]b = -l + 2[6]. (91) 

Since the only way for the b field to propagate is to decay 
into 5''^<I>", the scaling of the self energy of b is the same 
as the propagator of the composite operator ^"(J)". This 
observation results in 

[6] = 3- [*°$°]. (92) 

Thus rjh is known once we know the scaling dimension of 
the composite operator 

The scaling dimension of 5'°<I>° is obtained by adding 
the perturbation 

Ce~^ + (93) 

to the Lagrangian density and obtaining the flow of 5 by 
studying the Callan-Symanzik equation for the propaga- 
tor (*°(x)l>"(a;')). 

The Callan-Symanzik equation is essentially the "scale 
invariance" relation. 24 jv^gg^j. ^jjg fixed point, a change 
in scale has to be accompanied by the flow of some cou- 
pling constants and scaling of the flclds to maintain scale 




The result is 

(fc) = (^1+1 (^-| + ^ _ 1)^ log(|fc|/M) 

(99) 

from which we get 

A* = i + ^(^ + ^(e-i)). (100) 

Similarly, calculating the gauge invariant conduction 
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electron propagator 



by 



(101) 




(^fi"'(l + ^log(|fc|/M)) 



results in the gauge invariant 



(105) 



(102) 



The logarithmic infrared divergence of Eq. (101) can 
be interpreted as the scahng of the quasi-particle weight 
Z, defined by (^''^'') (k) = Z{if(:)-^ . The scahng form 
of Z which is consistent with Eq. (101) is 




(103) 



The vanishing of the Z factor as fc — > can in principle 
be seen in tunnelling experiments. This non-Fermi liquid 
behavior arises, because at the KI-ASL fixed point, there 
are very strong scattering of conduction electrons into the 
f-fermions. This is due to the fact that the b field which 
causes the scattering has become massless. 

The next step is to obtain ^^°c[)°^ (fc) in the prescnse 
of the perturbation 5^'^^'^ given by Eq. (93). This per- 
turbation is represented diagrammatically by 



(104) 



The diagrammatic expression for ^^I>°(()°^ is then given 



It results in 



-"l-^\og{\k\/n)]. (106) 



Callan-Symanzik equation [Eq. (97)] applied to the above 
propgator resutls in 



d5 _ 



The above relation together with Eq. (92) yields: 



Combine this with Eq. (91) to get 
1 / 16 



(107) 



(108) 



(109) 



Collecting what we have obtained for 7^ and r^f, given 
by Eqs. (86) and (109), together with the scahng relation 
given by Eq. (79), we arrive at the gauge invariant 



16 



37r2Ar' 



(110) 
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Note that the 1/N fluctuations have brought us further 
away from the generic 1/2 + 0(6) result one would obtain 
for a stable fixed point in an e expansion (which is an 
expansion around the gaussian action), assuming there 
exists such a stable fixed point. 



B. Staggered Spin Correlation 

We continue the study of the physical properties of our 
KI-ASL fixed point by obtainning staggered spin M ^ 
Sa — Sb correlation at the fixed point. The correlation 
is algebraic and our fixed point is an algebraic spin liquid 
fixed point. 

Since the decay exponent of {M{x) ■ M{x')) and 
{x)M~ {x')) is the same, it suffices to find the scaling 
dimension of M+. In terms of the field operators M^{x) 
is 



M+{x) = ^^{x)M<^^{x), 



(111) 



where is 4 x 4 matrix, which should be obtained from 
the relation between the field $ and the microscopic oper- 
ators developed in Sec. III. Physically 1 and 2 represents 
t and [. In 1/iV expansion they are 2 indices among 
many. M is obtained to be^^ 



M = l®az. 



(112) 



We should add this perturbation to the Lagrangian 
and study its flow. However $^(a;)A^$^(a;) mixes with 
^^(a;)A^\['^(x). In other words M+ is not a scaling op- 
erator. Having this in mind, we perturb the Lagrangian 
density by 

5Ce = U+ ^^{x)M^'^{x) + u+ ^^{x)M'^'^{x) (113) 

and study the Callan-Symanzik equation for the propa- 
gators (fc) and ($^$2^ (fc). The result (obtained 
in Appendix B) is 



d_ 

d£ 
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-2 



-2 



1 - 2 



(114) 



We emphasize that off-diagonal terms make 1 /N correc- 
tion to the scaling dimension of the staggered spin. This 
is because the scaling dimensions of \1/ and <& arc equal 
in the infinite N limit. This is similar to the importance 
of off-diagonal terms in the 1st order degenerate pertur- 
bation theory in quantum mechanics. 

The largest eigenvalue 1 -f- 2 (7 + ^73) /(Stt^TV) of 
above matrix dominates the scaling of Af+. The domi- 
nant scaling dimension of , denoted by [A/ +] , is then 
given by 



^A/^ 



2 (7 + 



31.1 
37r2Ar' 



The decay exponent of the correlation function results 
immediately: 



{M{x) • M(x')) oc 



1 



\X — X' 



■/\2[M-^ 



(116) 



(115) 



where the ellipsis denotes the decay with the larger ex- 
ponent coming from the other eigenvalue of the matrix 
in Eq. (114). 

Finally, we mention a comparison between our result 
and the Rantner-Wen studies. They studied spin cor- 
relations in their large-iV QED3 effective field theory in 
the Landau gauge. They were interested in the similar 
bilinear form as we have here with A4 = 1. In our cal- 
culation any bilinear form with the property [A4 , 7^] = 
has the same scaling dimension. To check our result we 
turned off all contributions to the exponent coming from 
terms involving The result [Af+]ASL = 2-32/(37r2Ar) 
obtained in Feynman gauge agrees with theirs obtained 
in Landau gauge. 2[Af+]ASL corresponds to the decay 
exponent of localized spin correlation when we are in 
the ASL phase, where the Kondo field is massive. As 
expected, we have [AT+Jasl < [-^~''], i-e. the Neel fluc- 
tuations are stronger in the ASL phase than at the KI 
transition, where the Kondo singlets arc forming. 

In the J < J'j^ regime the ASL is in the presence of 
the semi-metal (described by 5* field) and the gapped 
Kondo field b. The question may arise whether the semi- 
metal would destabilize the ASL. The semi- metal will 
not destabihze the ASL, because integrating out ^ and h 
fields in this phase will produce (retarded) four-fermion 
interaction term among $ fields. These four-fermion in- 
teractions are irrelevant at the ASL fixed point. 



VI. STABILITY OF THE FIXED POINT 

In this section we discuss the issue of the stability of 
the KI-ASL fixed point in the honeycomb lattice to 1st 
order in 1/A^. There are two separate issues one has to 
consider in the stabillity of this fixed point. 

In going to the continuum limit we have dropped the 
compact character of the gauge field. This is only justi- 
fied if the instanton events, where the flux of the gague 
field suddenly changes by 27r, do not proliferate. For 
Jk > Jk Kondo field b is condensed and we are 
in the Higgs phase. Gauge field is massive, its fluctua- 
tions are supperessed and instantons do not proliferate. 
Addressing the proliferation of instantons in the regime 
where the Kondo phase has disappeared is more subtle. 
We know from Polyakov's work that the compact pure 
gauge theory in 2+1 dimension is in the confining phase 
for all coupling constant, i.e. for all coupling constants 
the instantons proliferate. It has been argued recently 
that coupling the gauge field to massless Dirac fields re- 
sults in a deconfined phase for sufficiently large N }^ This 
is due to the fact that the massless matter field, in this 
case Dirac fields, suppress the fluctuations of the gauge 
field. We apply that argument here and assume that the 
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gauge field in the Jk < Jk regime is in a deconfined 
phase. At the KI-ASL fixed point, due to the presence 
of an extra massless matter field, i.e. the Kondo field b, 
the gauge field fluctuations are suppressed even more. 

Next we focus on the instanton-free sector. At the 
fixed point the gauge field a^, boson h and two fermion 
fileds 4*, all have scaling dimension 1 + 0{1/N). This 
can be seen by the form of propagators derived in Sec. IV. 
The large N scaling dimension of a^. b, 'J and (f> will 



limit our choices for relevant perturbations, 
such perturbations are discussed below: 



• b mass term: this term is already present in the 
microscopic bare Lagrangian. We have to tune Jk 
to make the effective mass of the b field to be zero. 
This perturbation is indeed relevant. 



and <J<i)°$°^ results in the following RG flow: 



d_ 



Ai 

A2 

(50 



= V 



■Ao- 

Ai 

A2 

6i 

S2 



(117) 



Different where the matrix anomalous dimension P is given by 



V 



1 



-74 32 32 



32 
32 
-2 

-4 
-4 



-74 32 
32 -74 

-4 -4 
-2 -4 
-4 -2 



-2 
-4 
-4 
-10 







-4 
-2 
-4 


-10 




-4" 

-4 

-2 





-10 
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• ^ and (f> mass terms: these mass terms either break 
the SU(N) flavor symmetry or the particle-hole and 
rotation symmetry. 

• b kinetic term of the form b* (9^ -I- ia^) b: this term 
violates the particle-hole symmetry. It has to be 
accompanied by 6(9^ — ia^j)b* with an equal co- 



efficient 
derivative 



But b*d^,b 



bd^b* 



dfj, {b*b) is a total 



2. 



b kinetic term of the form 6C = Ng\ [d^ + ia^) b\ 
this term is irrelevant by the power counting. That 
is 



dg 



0{1/N))i 



• Kinetic terms for 5* and <I> fields: these are the only 
potentially dangerous perturbations, because they 
are marginal in the N oo limit. Actually in the 
derivation of the Lagrangian density Eq. (61) we 
set Vc = Vf to make the action Lorcntz invariant. 
We have to consider perturbations that break this 
"tuned" Lorcntz symmetry and study their flow. 
By scaling argument the flow of these perturba- 
tions is marginal to first order in 1/N: marginally 
irrelevant, relevant, or exactly marginal. Our cal- 
culations results in the marginal irrelavance of these 
types of perturbations. We discuss this issue below. 

We break the Lorcntz symmetry by adding the pertur- 
bation 



to the Lagrangian density. The perturbation has to 
be considered both in and fields because of their 
mixing due to the exchange of the b field. Applying 
Callan-Symanzik equation to the propagators ^v]/"vfr"^ 



Eigenvalues of this matrix anomalous dimension are all 
negative except for one trivial zero eigenvalue. The zero 
eigenvalue is trivial since the corresponding perturbation 
can be absorbed by scaling "if and $ fields, shown in 
Appendix A. 

We knew from the mean-field result that there exist 
a fixed point in the parameter space. But our results 
goes further by ruling out a multi-critical fixed point. 
So based on our result, for large enough N, this quantum 
critical point will be achieved by tuning one parameter 
Jk. 



VII. CONCLUSIONS 

We have studied the Kondo-Heisenberg model on the 
honeycomb lattice at half filling. The main motivation 
for our study is that the Kondo band gap vanishes at 
a finite Kondo coupling constant (even as Jh — > 0) so 
that a novel kind of quantum phase transition becomes 
possible. 

In this paper we considered the transition of this 
Kondo insulator (KI) to a spin liquid phase which is later 
identified to be algebraic spin liquid (ASL). We developed 
a mean-field theory for such a transition using a fcrmionic 
representation for spins. To consider fluctuations about 
this mean-field systematically, we adopted a large N gen- 
eralization of our model by having N spin-fiavors rather 
than just {t,i}- We found that the KI-ASL transition 
is controlled by a stable Lorentz invariant fixed point. 

In the KI phase, the Kondo field b condenses and we 
are in the Higgs phase where the gauge field is massive. 
ASL phase is described by "decoupled" QED3 (coming 
from JnSi ■ Sj interaction) and the Dirac conduction 
electrons. They are "decoupled" in that the Kondo field 
that couples these two terms is now massive. The gauge 
field is deconfined in the ASL. 

We were able to calculate the exponent = I + 
16/(37r^A^) of the diverging length scale near the tran- 
sition. Interestingly, 1/A^ fluctuations have pushed the 
infinite N result further away from the result v = 1/2 
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one would get from the Landau-type expansion. We 
noted that at our fixed point the quasi-particle weight 
of the conduction electron vanishes which is indicative 
of non-Fermi liquid behavior. We also calculated the 
decay exponent of the staggered spin M = Sa — Sb 
correlation at the fixed point. This is not a scaling 
operator, since it mixes with m ~ sa — sb- The 
dominant decay exponent turned out to be 2[A/+] = 
4 - 4 (7 + \/73) /(Stt^TV) w 4 - 62.2/{3Tr^N). On the 
other hand, in the ASL phase M is a scaling operator 
with the decay exponent 2[Af+]ASL = 4 - 64/{3n^N). 
The jump in this exponent, as one gets to the KI tran- 
sition point, is due to the fact that there exists an ex- 
tra massless matter field at the transition (in this case 
Kondo field b), which modifies the exponent. This jump 
and its significance was discussed in a different context 
recently. 

The quantum phase transition between the KI and 
ASL does not break any physical symmetry. The KI- 
ASL transition is essentially a Higgs-deconfinement tran- 
sition. In other words, what changes at the transition 
point is the dynamics of the gauge field. This by itself is 
an interesting property of the transition we have studied, 
which is realized in a minimal Kondo-Heisenberg model. 

In connection to this work, we are currently pursu- 
ing the numerical study of this model; using a quantum 
Monte Carlo technique without the fcrmion sign problem. 
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APPENDIX A: RESTORATION OF THE 
LORENTZ INVARIANCE 



Here we will give more details regarding the irrelavance 
of the perturbations that break the Lorentz symmetry at 
our Lorentz invariant KI-ASL fixed point. It is tempting 
to break the Lorentz invariancc by perturbing only one 
field, or $. However "if and $ loose their own identity 
as b becomes massless and they mix together. Having 
this mixing in mind we consider the most general pertur- 
bation which is of the kinetic form for the ^ and $ fields, 
and is consistent with the gauge symmetry 



(Al) 

The corrections to fermion propagators and gauge vortex 
is represented diagrammatically: 




where we have used the notation 



= '7<t7<t- (A3) 



We evaluate the two propagators and ^<I>°<i>°) 

in the presence of the perturbation given by Eq. (Al). 



is given by 



Next we calculate ^^^^"^ which is given by 



where Sa-^^i, is given by 



5.;^. = 20J, + 2A, +4A^ + 4A,. (A6) 




(A4) 



where the sum over the space-time index cr, appearing 
in the above diagrams, is understood. Calculation of the 
diagrams — done in Fcynman gauge — results in: 



(fc) = -ir' ( 1 + iog(ifci//.) 



(A5) 






The above propagator has a contribution from the un- 
perturbed fixed-point action, which gives the gauge in- 
dependent scaling dimension for \1/ field 



A* = 1 



1 



(A7) 
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We obtain 

where Fg-;^^ is given by 



log(|fc|//i) 

log(|fc|/M) , 
157r2Ar ' 



(A9) 



-32A^ + 2(5^ + 45^ + 4,5^. (AlO) 



From Eq. (A9), we obtain the scaling dimension of the <I> 
field in Feynman gauge : 



1 



1 



(All) 



Then we use the Callan-Symanzik equation machinery 
by scaling the momentum k — > e~^k and applying 



| + (3-2A.) 
-| + (3-2A.) 



dC^ d 
dl dC„. 

dC„i d 
de dC„ 



(^«$") ^ 0, (A12) 
($'^1.'') = 0, (A13) 



where the set C includes the coupling constants that flow. 
For the perturbation under consideration it has 6 ele- 
ments: 

C = {Ao,Ai,A2,(5o,5i,(52}. 
Eqs. (A12), and (A13) will lead to the RG flow : 



d_ 



[Aol 




[Aol 


Ai 




Ai 


A2 

So 


= V 


A2 


Si 




Si 


.S2 




S2 



(A14) 



where the matrix anomalous dimension V is given by 



V 





"-74 


32 


32 


-2 


-4 


-4" 




32 


-74 


32 


-4 


-2 


-4 


1 


32 


32 


-74 


-4 


-4 


-2 




-2 


-4 


-4 


-10 










-4 


-2 


-4 





-10 







-4 


-4 


-2 








-10 



(A15) 



All the eigenvalues of this matrix anomalous dimension 
is nonpositive. There is one zero eigenvalue with the 
corresponding eigenvector 



A 
A 
A 
A 
A 
A 



The above shift can be absorbed by scaling 5" and <I> 
fields separately. The scaling is delicate, since in our RG 
scheme we have fixed the coefficient of (65'°$" + c.c.) to 
be 1. So the scaling of 5* and $ has to be accompanied 
by an appropriate scaling of b. But this - in general - 
will move us away from the critical point. For the above 
perturbation though, the scaling factor for 6 is 1: 



1 



1 + OiA^ 



So the above perturbation remains exactly marginal to 
all orders in 1/A^. 

It is interesting to note that the uniform shift 



is an eigenvector with a negative eigenvalue. Why should 
this perturbation flow? The absorbtion of the uniform 
perturbation in ^! and $ has to be accompanied by scal- 
ing the b field by 



^/TTA^/TTA 



1 - A + ©(A^ 



(A16) 



This scaling of the b field will move us away from the 
critical point. 



APPENDIX B: STAGGERED SPIN 
CORRELATION 

In this appendix we will find the algebraic decay ex- 
ponent of {M{x) ■ M{x')) at KI-ASL fixed point. We 
do this by obtaining the staggered spin M scaling di- 
mension. It suffices to find the scaling dimension of 
M+, since the decay exponent of {M{x) ■ M{x')) and 
{M^ {x)M~ [x')) is the same. The first step is to trans- 
late to the field theory language with the dictionary 
provided by Eq. 53. The expanded form is given by 



-ifA-. 
/!+ /b+ 



(Bl) 
(B2) 



Start with the expression for S\ — S\ + iS\ — f]^ f% 



which leads to 



S+{x) = ^'{x) 



(B3) 



The similar translation can be done for Sg. The final 
result is: 



and is then 



(B4) 



where A4 is 



M = 



/I \ 

1 

-1 

V -1/ 



^ 1(g) as. 



(B5) 



The matrix A4 has the following properties: 



Tr(7Wr) = 0, 



(B6) 
(B7) 



where F is a generic notation representing any number 
of 7^ matrices multiplied. The 1st property is useful 
since M passes through all the 7 matrices and the final 
result for the propagatros and (^E'^^'J'^) will have a 

simple matrix form ^Ad^t ^) . The second property 

causes any diagram, with matrix A4 left in a loop, to 
vanish. 

The scaling dimension of A/+ is obtained by looking 
at its relevant flow near the fixed point. Due to mixing 
we perturb the Lagrangian density by 



SCe = U+ ^\x)M^^ix) + u+ '^\x)M'^^ix), (B8) 



and study the Callan-Symanzik equation for the propa- 
gators ($^§2) and 

The vertices corresponding to these perturbations are 
represented by 



U+ M 



(B9) 



= u+ M 



The diagrammatic expression for propagators 
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The final answers in the Feynman gauge are as follows: 



J 



(fc) -r^A^r^ 



log(|fc|/M) 



C/+ 1- 



I 



log(|fc|/M) 



^^og{\k\/^,)\, 



TT^N 

2u+ 



log(|fc|//.) 



(B12) 
(B13) 



Then we use the Callan-Symanzik equation by scaling 
the momentum k e^^k and applying 



d_ 



-^ + (3-2A.) 





d 


d£ 


dCm. 


dCm 


d 


de 


dCra 



(B14) 



The set C includes the coupling constants that flow. For 
the perturbation under consideration it has two elements: 



The above Callan-Symanzik equations leads to the fol- 
lowing RG flow: 



d 






fl + 42_ 


dl 









-2 



371-2 AT. 



(B15) 



C^{U+,u+}. 
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